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By path integral Monte Carlo simulations we study the phase diagram of two - dimensional meso- 
scopic clusters formed by electrons in a semiconductor quantum dot or by indirect magnetoexcitons 
in double quantum dots. At zero (or sufficiently small) temperature, as quantum fluctuations of 
particles increase, two types of quantum disordering phenomena take place: first, at small values 
of quantum de Boer parameter q < 10~^ one can observe a transition from a completely ordered 
state to that in which different shells of the cluster, being internally ordered, are orientationally 
disordered relative to each other. At much greater strengths of quantum ffuctuations, at g ~ 0.1, 
the transition to a disordered (superfluid for the boson system) state takes place. 
PACS numbers: 

61.46+w: Clusters, nanoparticles, and nanocrystalline materials; 

68.65-|-g: Low-dimensional structures: structure, and nonelectronic properties; 

36.40Ei: Phase transitions in clusters 

The study of properties of mesoscopic systems is becoming very important due to the unabated reducing of the 
characteristic sizes of electronic devices. The progress in microlithography and semiconductor technology have enabled 
one to perform experiments with ultrasmall structures, containing as small as several electrons or excitons. Study of 
such systems has led to the development of one - electronics Q| and have provided essential progress toward a con- 
struction of different mesoscopic systems (see, e.g., [^) which can be considered as a base of nano- and optoelectronics 
of the future. 

In the present paper we consider the phase diagram of two-dimensional (2D) mesoscopic system of particles inter- 
acting through the Coulomb and dipole potential and confined by a quadratic potential. These models can be used 
when describing electrons in a semiconductor quantum dot (Coulomb cluster) or excitons in a system of vertically 
coupled quantum dots |^,^,^ (dipole cluster). Present experimental technique enables one to prepare the clusters of 
a given number of particles in both classical and quantum regimes ||^,|^ . Such possibility gives one an opportunity to 
explore a set of intriguing problems of the physics of mesoscopic systems. 

One of the most interesting features of mesoscopic classical clusters is the existence of the orientational disordering 
phenomenon (" orientational melting of small clusters" , see [|j - ) , when at temperatures many orders less than 
that of the total melting (of the disruption of cluster shells) mutual orientational disordering of different parts of a 
system takes place, pairs of shells rotate as a whole relative to each other loosing their relative orientational order. 

When the amplitudes of quantum fluctuations are comparable with the mean interparticle distance one can raise 
the question about the existence of a quantum analog of the classical orientational melting phenomena. One can 
expect that an increase in the strength of quantum fluctuations should lead to the lowering of orientational melting 
temperature Ts^s2 of ^ P^^ir of shells {si, S2} and, to all appearances, at some critical strength of quantum fluctuations 
"zero - point" orientational melting should take place. 

What defines the position of the point qs^s^ of zero-temperature orientational melting? How will be arranged the 
regions in which different pairs of a multishell cluster are orientationally ordered? What is the role of the statistics in 
this phenomena? In the present work we try to answer this questions with the help of the path integral Monte Carlo 
calculations. 

Mesoscopic clusters, properties of which we are studying, can be considered as two-dimensional system of N particles 
of mass m in a parabolic confinement potential of strength a. The general form of the corresponding Hamiltonian is: 

The forms of the interaction potential considered in this Letter are 1) the Coulomb interaction of electrons e in a 
semiconductor quantum dot U(rij) = jr^ and 2) the dipole interaction U{rij) = d^/r? . The latter interaction 
potential corresponds, e.g., to the system of indirect magnetoexcitons in a vertically coupled double quantum dot |^ 
or to electron system near metal gate that modifies interelectron interaction potential due to polarization The 
dipole moment d = he is defined by the distance h between the dots. 

Hamiltonian (|l|) can be reduced to a dimensionless form if all distances and energies are expressed in the units of 
ro and Eq — ar^, where tq — e^/'^/a^/'^ for a system of electrons (when the Coulomb cluster is considered; Cul); 
rg — d^/^/a^/^ for a cluster of magnetoexcitons (dipole cluster; D). Within the above units (|^) has the form: 
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^ = E (£i + + E + v^) + E f^N)' (2) 

i=l \ 2 / i<j 

U{n,)^{i/n,,Cul; 1/4, D} 

The strength of quantum fluctuations is controUed by the parameter q — H/ (m}/'^a^/'^r'^^ . The other dimensionless 
parameter which defines the state of the system is the dimensionless temperature T = khT/ (arg). 

In the path integral Monte Carlo approach (the multilevel block algorithm [|3| have been used) the properties of a 
given 2D quantum system {fi}, i = 1..N are estimated with the help of a fictitious 3D one {rf }, p = 0..P — 1 resulting 
from the discretisation of functional integrals. The required accuracy of this substitution, being controlled by the 
dimensionless parameter t = q/{PT), have been achieved by adjusting the number of levels of 3D system in such a 
way as to best satisfy the condition r = 0.3. 

The state of the cluster at each considered point {q; T} of dimensionless control parameters have been inferred 
from the analysis of the radial distribution function, radial mean-squared displacements of particles u^. Orientational 
disordering of a cluster have been tested with the help of the mutual orientational order parameter gsiS2 Ifl) . We 
have considered also mean - squared radial Z,- and angular fluctuations of trajectories in imaginary time, i.e. the 
measure of quantum " librations" of particles [p|Jl4| . 

Let us consider first the behavior of classical clusters (g — 0). The ground state configuration (at T <C 1) of a 
classical dipole system of 10 particles in a parabolic confinement has two distinct shells, the inner of which contains 3 
particles. The corresponding configuration can be designated as Z)io(3, 7). The insert of Fig. 1 shows that an increase 
in the temperature leads first, at T = T21 « 5.4 • 10~^, to the disappearance of the mutual orientational order of 
shells. Total disordering in the system, when radial shell order disappears, takes place at much more high temperature 
Tf w 0.01. 

It is obvious that in small clusters consisting of several shells a set of orientational melting phenomena will take 
place, each of them will appropriate to relative orientational disordering of different pairs of shells. This peculiarity 
of small clusters is well observed in Fig. 1 which represents temperature dependencies of radial fluctuations and 
of mutual orientational parameter gsiS2 of three - shell Coulomb cluster CUI25 the ground state configuration of 
which can be written as Cul25{3,9, 13). Note, that the data plotted in the figure reinforce the statement that the 
temperatures Ts-^S2 of the orientational "melting" of the pair of shells {si, S2} are functions of the particle distribution 
throughout shells (Ns-^jNg^) and are maximal when the shells involved are parts of ideal 2D hexagonal crystal (i.e. 
when (iV,,,iV,J = (3,9), (4, 10), (6, 12), ...). |Jl|] ^ 

As it follows from the Fig. 1, the temperature interval Ts-^S2 < T < Tf can be considered as a region in which 
shells {si, S2} of a cluster rotates relative to each other retaining their internal order and loosing the relative one. In 
2N - dimensional configurational space the system moves then along a narrow ravine on a potential energy surface, 
different points of this ravine correspond to different values of the mutual order parameter gsiS2- Fig- 1 shows that 
this movement is accompanied by a "breathing" of cluster shells to lead to sharp increase in the value of radial mean 
- squared fluctuations uf. at the points Tsj^..^ of orientational disordering. 

Let us consider how changes the state of clusters as the strength of quantum fluctuations is increased at constant 
temperature. Shown in Fig. 2a is the behavior of mutual orientational order parameter 1721 of dipole cluster Diq 
when the system is moved along line T = 3 • 10~^. A sharp variation in the value of (721 at the point 521 ~ 8 • 10~^ 
testifies about the quantum fluctuations induced transition from the oricntationally ordered (00, at q < 521) to the 
orientationally disordered but radially ordered state (RO, at g > (721). 

Analogous results take place for Coulomb cluster CUI25 and are plotted in the insert of Fig. 3a. The results shown 
in the figure correspond to the temperature T = 3 ■ 10~^ at which the orientational order of the first and the second 
shells of the cluster breaks down at (7 « 0.01. Quantum orientational melting of the second and the third shells takes 
place at sufficiently lower values of quantum parameter: ^32 « 8 • lO"'^. 

The dependencies of angular and radial Ir quantum librations of particles as functions of quantum parameter q 
is shown in Fig. 2b. In the insert of this figure a typical picture of imaginary - time trajectories (of their projections 
onto OXY plane) of the cluster is presented. The picture shows that quantum fluctuations are strongly anisotropic. 
An increase in the value of quantum parameter leads to a significant increase in the angular fluctuations at the point 
q2i of disappearance of an orientational order, while radial fluctuations Ir do not have any peculiarities in the region 
considered. 

Above results show that at the point of quantum orientational melting qsiS2 the characteristic scale of quantum 
motion becomes comparable with the mean angle between nearest particles of the pair of shells {si, S2}. A special 
feature of this type of quantum "melting" in mesoscopic systems is that at low temperatures main directions of 
particles fluctuations in corresponding classical system are defined by a narrow and high ravine formed by a potential 
energy surface. The presence of such ravine leads to a strong anisotropy of quantum fluctuations and sets off relative 
"quantum rotations" of a given pair of shells. It is also worth to note that orientational melting in both Coulomb and 
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dipole clusters take place in that region of dimensionless parameters {q; T}, in which the role of quantum statistics is 
unimportant (see the insert of Fig. 2b). 

Let us consider concisely how behave mesoscopic systems as the strength of quantum fluctuations is increased 
further, at g > 0.05. The simplest estimations show (see below) that in this region of phase diagram, as g > y/T, 
quantum exchanges of particles became to play an important role. In order to simplify our calculations, we have 
supposed the particles of the cluster to be bosons. Strictly speaking, in our case this can be applicable only for the 
dipole system of indirect magnetoexcitons in double quantum dots |^ . Electrons in a semiconductor dot obeys Fermi 
statistics and the model used collapses in the region of very strong quantum fluctuation. Meanwhile, the study of 
Coulomb boson system is of some methodical interest offering a clearer view of common features and regularities in 
behavior of boson clusters. 

In Fig. 3a the results of calculations of mean-squared radial fluctuations in dipole and Coulomb boson clusters 
at T = 5 • 10~^ are plotted vs. quantum parameter q. As it follows from the above analysis (see Fig. 1-2), in the 
considered region of control parameters all clusters are in orientationally disordered state. Therefore, a sharp increase 
in the value of testifies about the total "melting" in these systems. A comparison between the behavior of radial 
fluctuations (Fig. 3a) and that of the superfluid fraction (Fig. 3b) shows that this disordering corresponds to the 
transition from the radial ordered (RO) to the superfluid (SF) state. By way of illustration, in the insert of Fig. 3a 
we have plotted some typical pictures of imaginary - time trajectories of the system. Cyclic permutations of particles 
are well noticeable above the transition. As temperature Tf of the total disordering in appropriate classical systems 
is Tf « 0.01 and is approximately two times higher of the temperature to which the data in Fig. 3 belong, values gf 
of the total disordering can be taken as rather good estimations of locations of quantum phase transitions (at T — 0). 

By combining results considered above one can sketch the "phase diagram" of small quantum clusters on plane 
{q; r}, see Fig. 4. Of course, considering systems of such a small number of particles it is impossible to speak about 
the existence of distinct boundaries between states cited above, treating them as lines Tc{q) of phase transitions. But, 
as one can see from Fig. 1-3, an analysis of the system of even such a small number of particles enables one to reveal 
the presence of such regions as well as its arrangement in plane {q;T}. 

The mutual arrangement of different regions of the phase diagram is controlled by the following relations: 1) At small 
strength of quantum fluctuations the ratio of the temperature to the energy barriers for relative rotations of some pair 
of shells {S1S2} defines the region of orientational disordering T > Ts^s2 of this pair of shells. 2) Subsequent increase 
in the temperature (keeping q = const.) leads to the transition to the "classical liquid" (CL) state at T > Tf with the 
temperature Tf controlled by the ratio of the temperature of the system to the characteristic energy of interparticle 
interaction (i.e. by dimensionless temperature T in our units). 3) At T > T/, as the strength of quantum fluctuations 
is increased, and the thermal de Broigle wavelength becomes comparable with the mean interparticle distance in 
the liquid, the liquid / superfiuid transition takes place along line T{q) — q^. 4) Zero - temperature point g/ of a 
quantum melting, the transition of a cluster from a radial ordered to superfiuid (SF) state, is controlled by the ratio 
between the energy of quantum fiuctuations and the characteristic energy of interparticle interaction. 5) Comparing 
the energy of quantum fiuctuations of particles with the energy barrier for relative rotations of pair {si, S2} of cluster 
shells enables one to estimate the position of the quantum orientational melting point gsiS2- 

In conlusion, we note that the phenomena of quantum orientational melting can also take place near rare impurities 
of an infinite system, as in the vicinity of an impurity a shell structure can arise, resembling one of an isolated cluster. 

The work have been partially supported by Russian Foundation of Basic Researches, INTAS and by Sweden 
Academy of Science. 
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Figure 1 

Classical three - shell Coulomb cluster Cu/25, q = Q. The temperatures Ts^s2 of orientational melting of shells {si; S2] 
are marked by arrows. Shown also the temperature Tf of the total melting in the system, the transition from a radial 
ordered (RO) to a "classical liquid" (CL) state. 
Insert: Classical two - shell dipole cluster I?io, g = 0. 

Figure 2 

Quantum orientational melting of dipole cluster D^) at T — 3 • 10^^ and of Coulomb system CUI25 (in the insert of 
the Fig. 3a) at T = 5 • IQ-^. 

a) Mutual orientational order parameter 521 vs. quantum parameter q. 

b) A measure of quantum "librations" of particles: radial and angular variances of particles trajectories in imaginary 
time Iriq) and l,p{q). The region of abrupt changes in l^p coincides with that of the disappearance of orientational 
order. 

Insert: a typical picture of projections of imaginary - time trajectories onto OXY plane (the picture is appropriate to 
cluster Dio at the point {q; T} = {8 • lO^*; 3 • lO-^}). 

Figure. 3 
T = 5 • 10-3. 

a) Radial fluctuations vs. quantum parameter q. Inserts present spline - interpolations of OXY-projections of 
imaginary - time trajectories just before {q = 0.13) and after (q = 0.2) the 00/RO transition. 

b) Superfluid fraction i/g vs. q. 

Figure. 4 

The " phase diagram" of mesoscopic clusters (all notations are given in the text) . Dashed (solid) lines correspond to 
dipole cluster Diq (Coulomb system Cul2^). The region of transitions from orientationally ordered to orientationally 
disordered states is shown qualitatively on an enlarged scale. 
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